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Abstract 

The higher order susy partners of Schrodinger Hamiltonians can be exph- 

citly constructed by iterating a nonhnear difference algorithm coinciding with 

the Backlund superposition principle used in soliton theory. As an example, 

Qv I it is applied in the construction of new higher order susy partners of the free 

^D ■ particle potential, which can be used as a handy tool in soliton theory. 
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Q\ , Recent studies confirm that the higher order supersymnietric (susy) partners of 

^ \ Schrodinger Hamiltonians are most easily constructed by a simple algebraic tool named 

,-J^ \ intertwining technique [1]. One of the keys of this method is an algebraic nonlinear ex- 

pression which links solutions of different Riccati equations (see, e.g. [2-4]). In a previous 
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^ ■ paper [3] , we have studied the application of this method to the free particle potential. The 

'building blocks' of some of the resulting potentials are the well known soliton solutions of 

qh| the Korteweg-de Vries (KdV) equation: K^sech^[K(x — a)] and h? c'sc\<?[k{x — a)]. In this work 

we shall sketch the main steps of the approach in order to present some of the potentials 
derived in [3]. 

First, consider the intertwining relationship HiAi = AiHq, where the intertwiner Ai 
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^ is the first order differential operator Ai = -^ + (3i{x,e). All the available information 

concerning the Hamiltonians Hq = ~|;^ + K)(2;) and Hi = ~|^ + ^1(2^; e) is encoded in 
the beta function, which satisfies the Riccati equation 

-P[(x,e)+p'iix,e)=2[Voix)-e]. (1) 

The arbitrary integration constant e plays the role of a factorization energy. It is very simple 
to check that the potentials are related by the first order susy relationship 

Vi{x,e) = Vo{x)+f3[{x,e). (2) 

Equations (|I]) and (Q) are necessary and sufficient conditions for the Hamiltonians to be 
factorized as Hq — e = (1/2) AIAi, and Hi — e = {1/2) AiAi. Suppose now that Vq{x) is 
a known solvable potential with eigenvalues En and eigenfunctions t/'^, n = 0,1,2,... Let 
us assume that we have found a general solution of (|l]) for a given factorization energy 
ei 7^ En,Wn. Then, the potential V^i(x, ei) is also given [1-3]. The iteration of this procedure 



starts by considering now Vi(x, ei) as the known solvable potential and looking for a new 
one V2(x, ei,e) satisfying the second order susy relationship 

V2{x, ei, e) = Vi{x, ei) + P'2{x, ei, e). (3) 

Therefore, the new beta function must fulfill the Riccati equation 

- P'2{x,eue) + Piix,ei,e) = 2[Viix,ei) - e], (4) 

where e is again an arbitrary factorization energy. The corresponding solution is given by 

The finite difference expression (^ is a nonlinear superposition of two general solutions of 
(|I|), one for each factorization energy ei and e, transforming equation (|1|) into (|) by the 
change of Vo{x) by Vi{x, ei), and Pi{x, e) by /32(x, ei, e). This transformation can be used to 
link the higher order susy partners of Vo{x) with the first order superpotentials Pi{x, e), just 
by solving (^ for different values of the factorization energy e. For instance, providing n 
different general solutions of (|I]), one for each e^, k = 1, 2, ..., n, we are able to iterate n — 1 
times the algorithm (j^) acquiring a new beta function in each step, given by 

Pkix.ek) = -Pk-iix,ek-i)- -, ^^^ — -, r, k = 2,3,...n. (6) 

l3k-i{x, efc-i) - (3k-i{x, ek) 

We have adopted here an abreviated notation making explicit only the dependence of (3k on 
the factorization constant introduced in the very last step, keeping implicit the dependence 
on the previous factorization constants (henceforth, the same criterion will be used for any 
other symbol depending on k factorization energies). Therefore, given any initial potential 
Vb, the corresponding n-susy partner potential Vn can be writen as 

n 

Vn{x,en) = Vo{x) + Y.(3',,{x,ek), (7) 

fc=i 

provided that the master equations for (3k and Vk are given by 

- (3'^{x,ek)+ (3l{x,tk)=2[Vk-i{x,tk-i)-tk\, k = l,2,...,n, (8) 

Vk{x,ek) = Vk-i{x,ek-i)+ P'k{x,ek), k = 1,2, ...,n. (9) 

Now, let us stress that every general solution of the Riccati equation (|l]), for a given e, 
depends on an additional implicit integration parameter a, hence, the process acumulates 
as many of these integration parameters as many general solutions of (|lD have been used. 

Observe the coincidence of our nonlinear algorithm (|^) and the Wahlquist and Estabrook 
superposition principle expression (see equation (16) of [5]), derived from the Backlund trans- 
formation (BT) of the KdV equation Wt = Qwl — w^xx] subscripts t and x denote partial 
derivatives. The method has been typically used to generate new, multisoliton solutions 
Wi2, ...,w^n) of the KdV equation from a given one-soliton solution w = Wi oi the same 



equation. It is thus quite interesting that the vahdity of the same algorithm in the in- 
tertwining problem (supersymmetry) is much easier to demonstrate without worrying at 
all about the nonlinear equations! Moreover, its physical applicability in susy seems much 
wider. Thus, e.g., the singular solutions of KdV (singular water waves) would be of marginal 
physical interest. The singular potentials in the Schrodinger equation are not! Therefore, 
the possibility of reducing the n-th intertwining iteration to the multiple applications of the 
Backlund superposition principle means that n-susy could be a universal method generating 
the "multisoliton deformations" of any initial potential. 

We shall now focus on the vacuum case, presenting some simplifications which the method 
offers in deriving the n-susy partners for the potential Vq{x) = 0. In this case, the Riccati 
equation (|l|) has the general solution 

(3^{x,e) = -V2^cot[V2^{x-a)], (10) 

where a is an integration constant (in general complex). It is well known that the super- 
potential (p!0| ) gives four different first order susy partners of Vo{x) = by taking different 
values of e and a. This information is sumarized in Table I. 

TABLE I. The four different real superpotentials /3i comming out 
from din) , depending on the values of e and the integration parameter 
a. In each case S means singular, R regular, P periodic, and N null. 
The parameters a and b are arbitrary real numbers. 



Case e V2e a I3i{x,e) 



S e<0 i^y2\e\=iK a — k coth[K(2; — a)] 

27r 

R e<0 w2|e| = ZK —b — -— — k tanhlntx + b)] 

P e>0 \/2e = fc a -fc cot[fc(a; - a)] 

1 



N 



As an example, notice that the regular case (R) leads to the well known modified Poschl- 
Teller type susy partner V^{x,e) = — fi;^sech^[K(x + 6)], while the null case (N) leads to 
the potential barrier V^{x,0) = {x — a)^^. Now, in order to give an example of second 
order susy partner potentials V2{x,e), let us consider the superpotentials R and S as given 
in Table I. By introducing them in (^) and (^ we get 

V^ix €2) = -U^ - K^) ^1 ^^^^^ [^i(^ + ^)] + ^2 sech^ [/t2(x - a)] .^^. 

( —Ki coth [k,i{x + b)] + K,2 tanh [k2{x — a)] y 

The potential (|TTp has two finite wells which can be modulated by changing the values of 
Ki and K,2 under the condition k,2 < i^i- A Taylor expansion of (p!T]) shows a singularity at 
X = a when k,2 > i^i- The case K2 = ni gives a potential V2(x, ei) = 0. 



Let us remark that, for the periodic superpotentials Pi in Table I, equation (J^) leads to a 
natural classification of two kinds of potentials depending on the parity of n. For n even, the 
periodic superpotential /?i does not appear as a separate term in (0), affecting only one of 
denominators. The resulting susy partners have only a finite quantity of singularities. This 
fact has been used by Stalhofen [8] by constructing potentials with bound states embedded 
in the continuum. On the other hand, for n odd, the function jSi is a separate term in the 
sum ([^ and its global effect is not canceled by any similar term. The corresponding susy 
partners become singular periodic potentials. 

In conclusion, the nonlinear difference algorithm (^ allows the construction of higher 
order susy partners of any initial potential Vo(x), provided that a certain number of solutions 
of (P have been given. This finite difference algorithm generalizes the superposition principle 
reported in [5] extending its applications to the susy construction of new solvable potentials. 
In particular, the higher order susy partners Vn{x, e„) of the free particle potential represent 
a wide set of transparent wells in the terms discussed in [7-9] , as well as multisoliton solutions 
of the KdV equation as given in [5] . 
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